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CHAMJEL ESTIMATI0^7 IN A CELLULAR COMMUNICATION SYSTEM 



The present invention relates to a method and circuitry for 
channel estimation in a cellular communication system. 

The invention is parnicuiarly concerned v/ith a cellular 
communication system in which data is transmitted as a plurality 
of data symbols over a sequence of time slots. As is known, in 
a CDKiA system, data is encoded for transmission by modulating 
data symbols to be transmitted by a unique spreading code for 
each of a plurality of communication channels. Within each cell 
of a cellular communication system, spreading codes allow for a 
plurality of different mobile stations to communicate with a base 
station on selectively coded channels. 

When a signal is transmitted between a base station and a mobile 
station (either on the uplink or the downlink) , the signal 
receiving unit needs to establish from the signal which it has 
received some information about the communication path along 
which the signal has travelled. This is referred to herein as 
"channel estimation" and is-' carried out in a channel estimation 
unit which generates a channel impulse response. Various 
techniques are known for channel estimation. The channel impulse 
response is required in order to properly decode and demodulate 
incoming data. 

In an earlier proposed CDMA system explained in the Technical 
Report of lEICE (1996 - 02) entitled "Variable Rate Data 
Transmission on Single Code Channels in DS/CDMA" by Okomura et 
al, data is transmitted at variable transmission rates. Channel 
estimation is carried out from pilot symbols spaced periodically 
in the data stream. 

One difficulty associated with this proposal is that only the 
energy from the pilot symbols is available for channel 
estimation, and the signal to noise ratio of pilot symbols tends 
to be poor. Thus, filtering is required to get more reliable 
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estimates. In particular, in the case of mobile stations moving 
at high speeds, it may be difficult to apply a suitable filtering 
technique due to short channel coherence time. 

It is therefore advantageous to use. data symbols in addition to 
the pilot symbols for channel estimation. The present invention 
seeks to provide a technique for improving channel estimation 
using estimated data symbols. 

According to one aspect of the invention there is provided a 
method of channel estimation in a cellular communication system 
in which symbols are transmitted between a mobile station and a 
base station, the method comprising: 

receiving via a communication channel an incoming signal ; 
f^providing a first estimate for the channel impulse response 
of that communicatidn channel and utilising that first estimate 
to process the incoming signal to generate an estimated data 
symbol ; 

using said estimated data symbol to generate a soft output 
feedback decision; and 

using the soft output/ feedback decision to make a further 
estimate of the channel impulse response. 

The invention is particularly but not exclusively applicable to 
a CDMA system in which the incoming signal incorporates symbols 
spread by a spreading code. 

According to another aspect of the invention there is provided 
a method of channel estimation in a cellular communication system 
in which symbols are transmitted between a mobile station and a 
base station, the method comprising: 

receiving via a communication channel an incoming signal 
which incorporates symbols spread by a spreading code; 

providing a first estimate for the channel impulse response 
of that communication channel and utilising that first estimate 
to process the incoming signal to generate an estimated data 
symbol ; 
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making a further estimate of the channel impulse response 
by performing a correlation process between a conjugate version 
of the spreading code and the incoming signal, and multiplying 
the results of that process by a feedback decision based on said 
estimated data symbol as a weighting value for the further 
estimate of the channel impulse response. 

The incoming signal can include data symbols and pilot symbols, 
with a unique spreading code having been used respectively for 
piloc symbols and data symbols. 

According to the further aspect of the invention, the feedback 
decision may be a soft output or a hard output. 

A further aspect of the invencion provides circuitry for 
estimating a channel impulse response in a cellular communication 
system in which symbols are transmitted betv/een a mobile station 
and a base station, the circuitry comprising: 

a receiver for receiving via a communication channel an in 
coming signal which incorporates symbols spread by a spreading 
code ; y 

a channel impulse response generator for providing a first 
estimate for the channel impulse response of that communication 
channel ; 

a data symbol generator for generating an estimated data 
symbol using the first estimate of the channel impulse response; 

circuitry for providing a further estimate of the channel 
impulse response by performing a correlation process between a 
conjugate version of the spreading code and the incoming signal, 
and multiplying the results of that process by a feedback 
decision based on said estimated data symbol as a weighting value 
for the further estimate of the channel impulse response. 

A still further aspect of the present invention provides 
circuitry for estimating a channel in a cellular communication 
system in which symbols are transmitted between a mobile station 
and a base station, the circuitry comprising: 
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a receiver for receiving an incoming signal via a 
communication channel; 

a channel impulse response generator for providing a first 
estimate for the channel impulse response of that communication 
channel ; 

a data symbols generator for generating an estimated data 
symbol using that first estimate; 

a soft output feedback decision generator for using said 
estimated data symbol to generate a soft output feedback 
decision; 

wherein the soft output feedback decision is used to make 
a further estimate of the channel impulse response. 

The invention also contemplates a base station and/or a mobile 
station incorporating circuitry according to the first aspects 
of the invention. 

For a better understanding of the present invention and to show 
how the same may be carried into effect, reference will now be 
made by way of example to the accompanying drawings in which: 

Figure 1 is a block diagram illustrating schematically 
symbol estimation in a basic CDMA system; 

Figure 2 is a block diagram illustrating a cellular 
communication system; 

Figure 3 is a diagram illustrating a transmission time slot 
structure; 

Figure 4 is a block diagram illustrating channel estimation 
with feedback; 

Figure 5 is a block diagram of circuitry in a channel 
impulse response unit; 

Figure 6 is a vector space diagram; and 

Figures 7 and 8 are simulations of acquisition probability 
against signal energy per bit to noise ratio. 



Figure 1 is a block diagram pf a basic receiving system in a CDMA 
cellular communication network. Figure 2 is a block diagram 
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illustrating a context in which the present invention may be 
used. That is, a CDMA mobile communication system allows a 
plurality of mobile stations MSI, MS2 , MS3 to communication with 
a base station BTS in a common cell via RF signals on respective 
channels CHl, CH2, CH3 . These channels are discinguished from 
one another by the use of spreading codes in a manner which is 
known per se. 

Reverting now to Figure 1, signals incoming at an antenna 30 are 
received by an RF unit 28 and supplied to an analogue-to-digital 
(A/D) converter 32. The A/D converter takes samples of the 
incoming signals at intervals i to generate digital sampled 
values X. It will readily be understood that a signal may arrive 
at the receiving circuitry having experienced multi-paths with 
differing propagation delays. Each d;Lgital sampled value may 
include therefore components of a number of different symbols 
depending on the multipath effect. 

The digital signal is supplied from the A/D converter 32 to a 
channel impulse response unit 3 7 which provides an estimated 
channel impulse response W for use in processing the incoming 
signal. A reference code generator 22 supplies a reference code 
for despreading the incoming signal both to the channel impulse 
response unit 37 and to a rake receiver 36. The reference code 
is denoted C and is the unique spreading code for the 
communication channel on which the signal is received. Thus, it 
matches the code which was used to spread the signal prior to 
transmission at the base station or mobile station respectively. 
Although Figure 1 illustrates circuitry located at each mobile 
station, similar circuitry may also be utilised at a base 
station. 



The rake receiver 36 additionally ■ receives the input digital 
signal x and the estimated channel impulse response H. It 
generates an estimated transmitted signal u which is supplied to 
a decoder 31. The decoder estimates symbol values s from the 
estimated signal u which has been derived from the actual 
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incoming signal x, in a manner which is know per se. 

Operation of the rake receiver is known to a person skilled in 
the art but in any event is described in more detail in 
Application W096/24988 in the name of Nokia Mobile Phones 
Limited. 

The form in which data is transmitted is illustrated by way of 
example in Figure 3. The data symbols are transmitted in a 
sequence of time slots, each time slot including pilot symbols 
(PS) and coded data. Data is transmitted as a set of analogue 
symbols which are produced from originating digital data by i,q 
modulation. In a CDMA system the symbols are spread by a 
spreading code prior to transmission. In a typical system, each 
time slot has a duration of .625ms^. The pilot symbols are 
introduced at the beginning and end of each time slot in a 
coherent system. These symbols are readily recognisable and so 
the beginning and end of each time slot can be identified for 
synchronisation purposes. In the circuit of Figure l, the pilot 
symbols are additionally identified and used in the estimation 
of the channel impulse response by the channel impulse response 
37. 

Figure 4 illustrates an improved version of the circuitry of 
Figure l. In Figure 4, the estimated symbol s is fed back to the 
channel impulse response unit 3 7 via a feedback decision 
generator 20 which generates a feedback decision D from the 
estimated symbol s. That feedback decision D is used. to improve 
the estimate of channel impulse response generated by the CIR 
unit 37. An SIR block 29 receives the incoming signal samples 
X and the estimated signal u and calculates the signal to 
information ratio used to generate noise variance o- for the 
feedback decision generator 20. 

Figure 5 illustrates how the feedback decision D is- used to 
improve the estimated channel impulse response. Figure 5 
illustrates a CIR unit 37 which is capable of providing an 
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estimated channel impulse response H using both fedback estimated 
data symbols and incoming pilot symbols. The circuit comprises 
first and second matched filters 10,12 which are each used to 
apply the relevant despreading code to the incoming digital 
signal x^^. The reference code generator 22 generates three 
codes. A spreading long code LC is applied to both of the match 
filters 10,12. A spreading short code Wp for pilot symbols is 
applied only to the first match filter 10 and a spreading short 
code for data is applied only to the second match filter 12. 
The output signals 14,16 from the firsc and second match filters 
10,12 respectively are applied to an averaging circuit 18. The 
averaging circuit receives the decision feedback D and applies 
that to the signals 14,15 from the matched filters to generate 
the estimated impulse response H. 

The averaging circuit 18 includes a memory 40 for holding 
previous channel impulse responses and a weighting/averaging unit 
42 for taking an average of the earlier channel impulse responses 
and applying a weighting function according to the fedback 
decision value D, 

The correlation process which is used to establish the estimate 
of charmel impulse response can be expressed as .• 
H„ ( t ) = S D„»C* { i) x ( i+t ) =DnS C* ( i) x (i+t) (Equation 1) 
i i 

In the above equation, i is the sampling time index along the 
received signal, t is the resulting sampling time index of the 
correlation process, x(i+t) is the received data sample for that 
index, C* is the conjugate of the spreading code and is the n^^ 
decision feedback where n is the time index of the impulse 
response. It will be apparent that the decision feedback D may 
be in respect of pilot symbols and/or data symbols. 

In use, an initial estimate for the channel impulse response H 
would be established by the CIR unit 3 7 on the basis only of 
known incoming pilot symbols in a particular time slot in a 
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manner which is known per se. Alternatively, an initial estimate 
can be made using Equation (l) without the weighting function D. 
Once a preliminary estimate of the channel impulse response has 
been established, this is held in the memory 40 and supplied to 
the rake receiver 35 for subsequent signal processing of data 
symbols in that time slot. The first data symbol to be estimated 
is then subject to feedback to provide a feedback decision D so 
that a new estimate for the channel impulse response K can be 
determined using the procedure outlined above and as illustrated 
in Figure 5. The spreading code conjugate C* in equation 1 is 
derived from merging the conjugates of and LC (for data} or W 
and LC (for pilot symbols) . Then, an average is taken of the new 
estimate and the preliminary estimate. Thus, at each estimated 
symbol an improved channel impulse response is obtained for 
subsequent processing. The average can be taken coherently 
within any coherent period, that is part of a slot, one whole 
slot or a multiple of slots. Alternatively, a non-coherent 
average can be taken over any period. 

It will be apparent from Equation 1 that the correlation process 
can effectively be implemented without decision feedback as if 
all the transmitted data symbols are ones. The transmission 
feedback D acts as an average weighting factor to determine the 
estimated channel impulse response H. 

Another aspect of the implementation discussed herein is to 
determine the nature of the decision feedback D generated by the 
feedback decision generator 20. Figure 6 illustrates the theory 
underlying one proposal which is to utilise a soft symbol 
decision feedback (SDPB) process by generating a soft feedback 
decision value D. Figure 6 illustrates a vector diagram of 
signal space. Vector A is a transmitted data signal. Vector A' 
is the received data. Vector A is the actual value we expect 
from the transmitted data signal since the rake receiver cannot 
be perfect. The received data symbol can be written as: 



A' =a+n 



(Equation 2) 
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where n is the noise vector. 
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it is natural to assume that both the signal and the noise are 
zero mean. We have: 

E{a} =0 
E{n} =0 

E{A'}=E{a+n} =0 (Equation 3) 

In a binary (+1/-1) environment, we have also: 
E{AA*} =1 

E{A'A'*} =E{ (a+n) (a+n) *) 
=E{a a*}+ E{nn*} 

=E{a a*}+ cr^ ^ (Equation 4) 

where is the noise variance. 

With these statistics of the signal, the soft output can be 
adjusted to satisfy the conditions for feedback if the noise 
variance is known beforehand. For a simple solution, we suppose 
that the signal loss is minimum and can be neglected: 

|a-A| -> 0 

Then equation (3) becomes: 

E{A'A'*} =E{a a*}+ 

=E{A A*}-r a' 

= 1 + (Equation 5) 

Using this analysis, the feedback decision generator can generate 
a soft value D satisfying the following requirements: 

E{D} = 0 (Equation 6) 



E{DD*}- = 1 + 



(Equation 7) 
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As an alternative, a hard decision feedback may be utilised by 
performing a binary cut-off, that is to provide either +i or -i 
to the input of the channel impulse response unit for each 
estimated symbol S. 

The effectiveness of both solutions is illustrated in Figures 7 
and 8. In figures 7 and 8, the x axis represents signal energy 
per bit to noise ratio {E^/NJ and the y axis represents the 
scale of acquisition probability - that is the likelihood that 
all signal propagation paths have been captured in the received 
signal as estimated using the C.I.R. Figure 7 is the simulation 
results with the best possible channel estimation, and Figure 8 
is the simulation results in which the channel information of a 
slot is obtained only from the pilot symbols in the same slot. 
These simulations have been generated^based on the MATLAB program 
included as Annexe 1. The simulation results compare the 
acquisition probabilities of the impulse response of different 
approaches of the decision feedback compared with no decision 
feedback. In both cases, it can be seen that the best results 
are given by a soft decision feedback for both pilot and data 
symbols. This is the top/ line in both simulations. The next 
line is for soft feedback decisions relating only to data 
symbols, that is without pilot symbol feedback being taken into 
account . 

The next line is for hard data feedback decisions in respect of 
data symbols only. 
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The lowest line is where there is no decision feedback. 
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ANNEXE 1 

MATLAB program 

% Soft Decision Feedback 

% ARIB specification 

% QZD/NTC - FEB. 98 - Oulu 

clear, rand (' seed' , o ) , randn ( ' seed' , 0 ) 

Dat_Weight = 0.5; SNR = -30:3:0; Fingers = 10; 

h= [0 . 7+j *0.1,0,0,0,-O.S-j*0.3,0,0,0,0,0.3+j*0.5,0,0, 

0.1-j*0.13,0,0,0,0,0,0,G.3-j*0.1,0,G,0,0.5+j*0.5! 

0,0,Q,-0.5-j*0.3,0,Q,0,0,0.3+j*0.5,G,0,0.1-j*0.13 

0,O,O,0,0,0,0.3-j*0.l3; 
h=h/sqrt (h*h' ) ; 

AcqError = zeros (size (SKR, 2 ) , 4); 

for kkk = 1 :size (SNR, 2) 

snr = SNR (kkk) % SNR in dB 

for nnn=l;10 0 0 
DAT = (rand(l,40) >0.5) *2-l; 
% CODE 

CL_I= (rand (1,25 SO) >0.5) *2-l; CL_Q= { rand ( 1 , 25S0) >0.5) *2-l; 
C = CL_I + j*CL_Q; 

w=l;for k=l:6;w=[w W;W - w] ; end ; wO =w (l , : ) ; wl =w ( 2 , : ) ; clear W; 
WO=r];Wl=[]; for k=l:40, W0= [WO, wO] ;W1= [Wl, wi] -end 
CO = wO .*C; CI = Wl .*C; 
% Spreading 

D=[]; for k=l:40, D- [D,DAT{k) *wl] ;end % Data 
P=ones (size (D) ) ; % pilot 

XI = (D + j*P) .* C; 
% Channel 

x2 = conv(xl,h) /2; 
% Noise 

nA=l/lO" ( snr/ 20);fork=l:size(x2,2);n(k) =nA*randn*exp ( j *rand*2 *pi ) ; end 
X = x2 + n; ' 

%%%%%%%%%%%%%%%%%%% MF %%%%%%-%%%%%%•%%%■%%%% 
IRPilot = zeros (39, 64); IRData = zeros(39, 6^); 
fork=l:39, l=(k-l)*64; 
for m = 1:64 

IRPilot (k,m) = IRPilot(k,m) + x ( l+ra- 1+ (1 : S4 ) ) *C0 ( 1+ (l : 64 ) ) ' 
IRData(k,m) = IRData(k,m) + x ( l+m-l+ ( 1 : 64 ) ) *C1 ( 1+ (1 : 64 ) ) ' ; 
end 

end 

IRPilot = -j*IRPilot/64; 
IRPilotHard=sum{IRPilot (1:28, : ) ) /28; 
IRData = IRData/64; 

%%%•%%%■%%%%%%%%%%%%% H' %%%%%%%%%%%%•%%%%%%%%%%%%%% 
a = IRPilotHard .* conj ( IRPilotHard) ; 
H = zeros (size (h) ) ; 
for k = 1: Fingers, findf = 0; 
while ( findf < 1) 

[n, m] =max (a) ; a(ra)=0; 

if (m< = size (h, 2) ) , H (m) =IRPilotHard (m) ; findf=2; end 

end 

end 

%%%%%%%%%%%%%%%%%%% decision %%%■%%%%%•%%%%%%%%%%% 

hl=find(abs (h) >0 . 01) ; 

h2=hl; % Idea 

%h2=f ind(abs (H) >0 . 01) ; % blind 

PilotSoft = zeros (1,28); DataSoft = zeros (1,40); 
for n=l: Fingers 

a = x(h2 (n) :size (CO, 2)+h2 (n) -1) .* conj (CO) * conj (H (h2 (n) ) ) ; 
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for k=l:28 

PiloCSoft(k> = PilotSoft(k) + imag (sum (a ( (k- 1 ) *64 + 1 : k*64 ) ) ) ; 
end 

a = x(h2 (n) :size(Cl,2)+h2(n) -1) .*conj(Cl) * conj (H {h2 (n) ) ) 
for k=l:40 

DataSofC(k) = DataSoft(k) + real ( sum (a ( (k-1 ) * S4 + 1 : k*64 ) ) ) ; 
end 

end 

PiloCSoft = PilotSof t/64; 
DataSoft = DataSof t/64 
DataHard = sign (DataSof t) ; 

%%%%%%%%%■%%%%■%%%%%% Averaging %%%%%%%%%%%%%%%%%%% 
a = 0 ; b = 0 ; c = 0 ; 

IRHard=zeros (l, 64) ,- IRSoft = zeros (1, 54) ; IRSof CD=zeros (1, 64) ; 
IRSof ter=zeros (1, 64) ; IRSof tero=zeros ( 1 , 64) ; 
for k=l:39 

IRHard=.IRHard+IRData (k, : ) *DataHard (k) ; a=a-rabs (DataHard (k) ) 
IRSof t=IRSoft+IRData (k, : ) *DataSof t (k) ; b^b+abs (DataSoft (k) ) ; 

end 

for k=l:28 

IRSof tp=IRSoftp+IRPilot (k, : } *PilotSof t (k) ; c=c+abs 
(PilotSof t (k) ) ; 

end 

IRSoftp = IRSoft*Dat_Weight/b + IRSoftp * ( l-Dat_Weight) /c ; 

IRHard = IRHard*Dat_Weight/a + IRPilotHard* (l-DaC_Weight) ; 

IRSoft = IRSoft*Dat_Weight/b + IRPilotHard* ( l-Dat_Weight) ; 

%************ acq probability ************ 

a = IRPilotHard .* con j (IRPilotHard) ; 

for k=l : size (hi, 2) , [n , m] =max (a) ; h2 (k) =m; a (m) =0 ; end 

c=0; for k=l :size (hi , 2) , f or n=l : size (hi , 2 ) 

if (h2 (k) ==hl(n) ) , c=c+l; h2(k)=0; end 
end, end 
AErrNoDFB = size (h.1 , 2 ) -c ; 
a = IRHard .* con j ( IRHard) ; 

for k=l: size (hi, 2) , [n,m] =max{a) ;h2 (k) =m;a (m) =0;end 
G=0; for k=l : size (hi , 2) , f or n=l : size (hi , 2 ) 

if (h2 (k) ==hl (n) ) , c=c + l; h2(k)=0; end 
end, end 
ASrrHardDFB = size (hi, 2) -c; 
a = IRSoft .* conj (IRSoft) ; 

for k=l :size (hi, 2) , [n,m] =max(a) ;h2 (k) =m;a(rn) =0;end 
c=0; for k=l :size (hi, 2) , for n=l size (hi , 2 ) 

if (hi (k) ==hl (n) ) , c=c+l; h2{k)=0; end 
end, end 
AErrSoftDFB = size (hi, 2) -c; 
a = IRSoftp .* conj (IRSof tp) ; 

for k=l :size (hi, 2) , [n,m] =max(a) ;h2 (k) =m;a (rn) =0 ;end 
c=0; for k=l :size (hi, 2) , for n=l : size (hi , 2 ) 

if (h2(k)==hl(n)) , c=c+l; h2(k)=0; end 
end, end 
AErrSoftPDFB = size (hi , 2 ) -c ; 
AcqError (kkk, : ) =AcqError (kkk, :) + ... 

[AErrNoDFB, AErrHardDFB, AErrSoftDFB, AErrSoftPDFB] 
end % nnn 
end % kkk 
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